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ABSTRACT 



Aims. We study turbulent transport coefficients that describe the evolution of large-scale magnetic fields in turbulent convection. 
Methods. We use the test field method, together with three-dimensional numerical simulations of turbulent convection with shear 
and rotation, to compute turbulent transport coefficients describing the evolution of large-scale magnetic fields in mean-field theory 
in the kinematic regime. We employ one-dimensional mean-field models with the derived turbulent transport coefficients to examine 
whether they give results that are compatible with direct simulations. 

Results. The results for the a-effect as a function of rotation rate are consistent with earlier numerical studies, i.e. increasing magni- 
tude as rotation increases and approximately cos 6 latitude profile for moderate rotation. Turbulent diffusivity, r]t, is proportional to the 
square of the turbulent vertical velocity in all cases. Whereas rjt decreases approximately inversely proportional to the wavenumber 
of the field, the a-effect and turbulent pumping show a more complex behaviour with partial or full sign changes and the magnitude 
staying roughly constant. In the presence of shear and no rotation, a weak a-effect is induced which does not seem to show any 
consistent trend as a function of shear rate. Provided that the shear is large enough, this small a-effect is able to excite a dynamo 
in the mean-field model. The coefficient responsible for driving the shear-current effect shows several sign changes as a function of 
depth but is also able to contribute to dynamo action in the mean-field model. The growth rates in these cases are, however, well 
below those in direct simulations, suggesting that an incoherent a-shear dynamo may also act in the simulations. If both rotation and 
shear are present, the a-effect is more pronounced. At the same time, the combination of the shear-current and fi x ./-effects is also 
stronger than in the case of shear alone, but subdominant to the a-shear dynamo. The results of direct simulations are consistent with 
mean-field models where all of these effects are taken into account without the need to invoke incoherent effects. 
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1. Introduction 

The solar magnetic field is thought to arise from a complicated 
interplay of turbulence, rotation, and large-scale shear flows 
(e.g. Ossendrijver 2003 and references therein). Whilst numer- 
ical simulations of simple systems using fully periodic boxes 
and externally forced idealised flows exhibiting large-scale dy- 
namos have been around for some time (e.g. Brandenburg 
I200T1 I2005al Brandenbu rg et al. I200T1 Minin ni et a l . [2005l 
Brandenburg & Kapyla [20071 Yousef et al. |2008al |2008bl 
Kapyla & Brandenburg 12009b and dynamos driven by the 
magnetorotational instability exhibit large-scale dynamos (e.g. 
Brandenburg et al. 1995, Hawley et al. 1 19961 1, convection sim- 
ulations have not been able to produce appreciable large-scale 
magnetic fields until recently (Rotvig & Jones 2002 ; Browning 
et al. 120061 Brown et al. I2UU71 K apyla et al. |2008l hereafter 
Paper I; Hughes & Proctor 2009). The main ingredient miss- 
ing in many earlier simulations was a large-scale shear flow and 
boundary conditions which allow magnetic helicity fluxes out of 
the system. Indeed, the shear flow plays a dual role in dynamos: 
it not only generates new magnetic fields by stretching, but it 
also drives magnetic helicity fluxes along constant isocontours of 
shear which can allow efficient dynamo action (Vishniac & Cho 
2001, Brandenburg & Subramanian 2005; Paper I). Recently, 
however, large-scale dynamos have also been found from rigidly 
rotating convection simulations without shear (Kapyla et al. 
20093l. 



Although large-scale magnetic fields can clearly be obtained 
from simulations, the origin of these fields in many cases (e.g. 
Yousef et al. l2008all2008bl Paper I; Hughes & Proctor|2009) is 
still uncertain. In the mean-field framework (e.g. Moffatt 1978; 
Parke r [19791 Krause & Radler [19801 Riidiger & Hollerbach 
2004), the dynamo process is described by turbulent transport 
coefficients that govern the evolution of large-scale magnetic 
field. The evolution equation for the large-scale part is obtained 
from the standard induction equation by decomposing magnetic 
and velocity fields into their mean and fluctuating parts, i.e. 
B = B + b,U = TJ + u, which leads to 



Q-Q 

— = Vx(UxB + £- 77/zo J) 



(1) 



where rj is the molecular magnetic diffusivity, J = /i V x B 
is the current density, and //o is the vacuum permeability. The re- 
maining term, £ = u x b, is the electromotive force describing 
the effects of small-scale turbulence on the evolution of mean 
fields and can be represented in terms of the mean fields and 
their derivatives 



£ i = a. 
where a 



(2) 
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ij and rjijk are tensorial coefficients, commas denote 
partial derivatives, and summation over repeated indices is as- 
sumed. Expression (|2]i is valid if the mean fields vary slowly in 
space and time. 

Whilst mean-field models have been quite successful 
in reproducing many aspects of the solar magnetism (e.g. 
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Ossendrijver 2003 ), they have often been hampered by the poor 
knowledge of the turbulent transport coefficients which could 
only be computed analytically using unrealistic or unjustified 
approximations, such as first order smoothing (FOSA). More re- 
cently, numerical models of convection in local Cartesian geom- 
etry have been employed to compute some of these coefficients 
in more realistic setups (Brandenburg et al. 1990; Ossendrijver 
et al. I200T1 [20021 Giesec ke et al. I2U051 Kapy la et a l. I2006al 
Cattaneo & Hughes [20061 Hughes & Cattaneo l2008l i. To date, 
however, only coefficients relevant for the ay term in Eq. (O 
have been determined from convection simulations. This is due 
to the limitations of the method used where a uniform magnetic 
field is imposed and the resulting electromotive force is mea- 
sured. Furthermore, if the Lorentz force is retained in the simu- 
lations, dynamo-generated magnetic fields may grow to satura- 
tion, leading to quenching even if the imposed field is weak. At 
large magnetic Reynolds numbers such quenching can be very 
strong if there are no magnetic helicity fluxes, suggesting there- 
fore small values of a even for weak imposed fields. 

During recent years an improved scheme of extracting tur- 
bulent transport coefficients has appeared which is referred to 
as the test field method (Schrinner et al. 120051 [2007b . In the test 
field method the velocity field of the simulation is used in a num- 
ber of induction equations, which all correspond to a given set 
of large-scale test fields which do neither evolve nor react back 
onto the velocity field. The test fields are orthogonal so the co- 
efficients can be obtained by matrix inversion. This method has 
been used successfully in setups where the turbulence is due to 
isotropic forcing without shear (Sur et al. 120081 Brandenburg et 
al. [2008 h) and with shear (Brandenburg et al. |2008at Mitra et 
al. 12009), respectively. Moreover, the method has been used to 
extract dynamo coefficients from more realistic setups where the 
turbulence is driven by supernovae (Gressel et al. 2008 ) and the 
magnetorotational instability (Brandenburg 2005b, 2008 ). 

In the present paper we apply the method for the first time 
to convection simulations. We also seek to understand the dy- 
namos reported in Paper I by applying the derived coefficients in 
a one-dimensional mean-field model. In the case of convection 
with rigid rotation it is likely that the large-scale fields are due 
to the turbulent a-effect that is present in helical flows (Kapyla 
et al. |2009al >. However, when shear is present, there are vari- 
ous mechanisms that can generate large-scale fields: in helical 
flows a finite a-effect (e.g. Radler et al. 2003 1 Radler & Stepanov 
2006; Rudiger & Kitchatinov 2006 ) with shear can excite a clas- 
sical afl or a-shear-dynamo (e.g. Brandenburg & Kapyla 120071 
Kapyla & Brandenburg 2009). Even if the mean value of a is 
zero, strong enough fluctuations about zero in combination with 
shear can drive an incoherent a-shear dynamo (e.g. Vishniac 
& Brandenburg 1 1 9971 Proctor 2007 ). Finally, the shear-current 
(Rogachevskii & Kleeorin 2003 2004; Kleeorin & Rogachevskii 
2008 ) and £1 x J (Radler [19691 Radler et al. |2003l Pipin |2008l 
effects may operate even in nonhelical turbulence. If both rota- 
tion and shear are present in the system it is not obvious how 
to distinguish between the shear-current and fi x J effects. In 
the present paper we are able to extract the relevant turbulent 
transport coefficients responsible for most of these processes and 
determine which one of them is dominant in the different cases 
with the help of a one-dimensional mean-field model. In order 
to facilitate comparisons between the mean-field models and the 
direct simulations presented in Paper I, we use identical setups 
and overlapping parameter regimes as those used in Paper I in 
the determination of the transport coefficients. 



2. Model and methods 

The setup is similar to that used by, e.g., Brandenburg et al. 
([19961 Ossendrijver et al. J2001II2002I ). and Kapyla et al. d2004l 
2006a) and in Paper I. A small rectangular portion of a star is 
modelled by a box situated at colatitude 8. The coordinate sys- 
tem is such that (x, y, z) corresponds to (9, <fi, r) in a spherical 
coordinate system. The dimensions of the domain are in most 
cases (L x , L y , L z ) = (4, 4, 2)d, where d is the depth of the con- 
vectively unstable layer, and it is also used as our unit length. 
The box is divided into three layers, an upper cooling layer, a 
convectively unstable layer, and a stable overshoot layer (see be- 
low). The following set of equations for compressible hydrody- 
namics is being solved: 



V\np 
Vt 



-V-C7, 



(3) 



^- = -SU x y--Vp + g-2ftxU + -V -2upS, (4) 
Vt p p 



^ = - P -V-U + l -V-KVT + 2^- C ^, 

Vt p p t(z) 



(5) 



where V/Vt = d/dt + (U + U {S) ) • V, and U {S) = (0,Sx,0) 
is the imposed large-scale shear flow. The kinematic viscosity is 
given by v, p is the density, U is the velocity, and g = —gzis the 
gravitational acceleration. The fluid obeys an ideal gas law p = 
pe(7 — 1), where p and e are the pressure and internal energy, 
respectively, and 7 ~ cp/cy = 5/3 is the ratio of specific heats 
in constant pressure and volume. The internal energy is related 
to the temperature via e = cyT, and K is the heat conductivity. 
The rate of strain tensor S is given by 



■ u 



(6) 



The last term of Eq. (|5]l describes cooling at the top of the do- 
main, where t(z) is a cooling time which has a profile smoothly 
connecting the upper cooling layer and the convectively unstable 
layer below. 

The coordinates {z\, Z2, 23, Z4) = (—0.85, 0, 1, 1.15)rf give 
the vertical positions of the bottom of the box, the bottom and 
top of the convectively unstable layer, and the top of the box, 
respectively. We use a K(z) profile such that the associated hy- 
drostatic reference solution is piecewise polytropic with indices 
(mi, 77i2, "73) = (3, 1, 1). The cooling layer near the top makes 
that layer nearly isothermal and hence stably stratified. The bot- 
tom layer is also stably stratified, and the middle layer is convec- 
tively unstable. 

Stress-free boundary conditions are used for the velocity, 



U x 



u z = o. 



(7) 



In the absence of shear the x and y directions are periodic 
whereas if shear is present, shearing-periodic conditions are used 
in the x direction. A constant temperature gradient is main- 
tained at the bottom of the box which leads to a steady influx of 
heat due to the constant heat conductivity. The simulations were 
made with the PENCIL CODlfl which uses sixth-order explicit 
finite differences in space and third order accurate time stepping 
method. Resolutions of up to 256 3 mesh points were used. 



http : / /www . no rdit a . org/ software/pencil-code/ 
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2.1. Units, nondimensional quantities, and parameters 

Dimensionless quantities are obtained by setting 

d = g = p = c P = /i = 1 , (8) 

where po is the density at z 2 . The units of length, time, velocity, 
density, entropy, and magnetic field are then 



[x]=d, [t] = Vd/g, [U] 

[s] = c P , [B] 



dg , [p] = Pa 



gpopo 



(9) 



The simulations are then governed by the dimensionless num- 
bers 



Pr = — 

Xo 



Re 



vkr 



Ra 



g^_ 



_L_ds 

cp dz 



where \o — K/(p m cp) is the thermal diffusivity, kf = 2w/d 
is an estimate of the wavenumber of the energy-carrying ed- 
dies, and p m is the density in the middle of the unstable layer at 
z m = \{z3—Z2). Our choice of fcf is somewhat arbitrary because 
it is difficult to define a single length scale which would describe 
the flow in a highly inhomogeneous system such as stratified 
convection. The vertical extent of convective cells, however, is 
almost always of the order of the depth of the convectively un- 
stable layer which suggests that d could be used as the length 
scale describing convection. In the nonrotating case, this is also 
close to the horizontal size of the convective eddies. The entropy 
gradient, measured at z m in the non-convecting initial state, is 
given by 



J_ds\ = V- V ad 
cpdz-/ Hp 



(11) 



with V ad = 1 - I/7 and V = (<91nT '/ 'dlnp) Zm , and H P is the 
pressure scale height. 

The amount of stratification is determined by the parameter 



Co 



(7 - l)eo 

gd 



(12) 



where eo is the internal energy at Z4. We use £0 = 1/3 in all 
models. 



2.2. The test field method 



We employ the test field method (Schrinner et al. 120051 120071) . 
which is implemented into the PENCIL Code, to determine tur- 
bulent transport coefficients. The uncurled induction equation in 
the shearing box approximation can be written in terms of the 
vector potential in the Weyl gauge as 



DA 

= -SA y x + U x B - i]p J, 



(13) 



where D/Dt = d/dt + Sxd/dy, A is the magnetic vector po- 
tential, and B — V x A is the magnetic field. The relative im- 
portance of magnetic diffusion over viscous and inertial forces 
can be characterized respectively in terms of magnetic Prandtl 
and Reynolds numbers 



v 

Pm = — . 



Rm 



PmRe. 



(14) 



77 rjkf 

In most cases we use Pm = 5 and Rm w 35, see Table [TJ 
When we vary Rm in the range from roughly 1 .5 to 150, we keep 



Re w 15 and vary Pm in the range 0.1 — 10. We decompose the 
fields into their mean and fluctuating parts according to 



U = U 



B = B + b, J = J 



(15) 



where the overbars denote a horizontal average and lowercase 
quantities denote fluctuations around these averages. The equa- 
tion for the mean vector potential is then 



DA 
~Di 



~ S AyX 



UxB+uxb— r)(ioJ. 



(16) 



Subtracting ( TT~6b from (13[ gives an equation for the fluctuating 
field which reads 

Da 



(10) ~Dt 



-Sa,yX + U xb+uxB+uxb—ux b — r]fi j. (17) 



Instead of using the actual mean fields B in this equations, 
they are replaced by orthogonal test fields B ' and a separate 
Eq. (l% is solved for each one of them. Here we follow the 
same procedure as in Brandenburg et al. (2008a) and Mitra et 
al. (J2009) and limit the study to mean magnetic fields that de- 
pend on z only. We use test fields 



B 



£?o(cos kz, 0, 0), B = _B (0,cosA:z,0), 



B ls = Bo(sin/fcz,0,0), B 2s = B (Q, sin kz, 0), 



(18) 
(19) 



where k is the wavenumber of the test field. In most models we 
use kjk\ = 1, where k\ = 2ir/L z . The electromotive force can 
be written as 



j • 



(20) 



where r/u = 77^23 and 7^2 = —77^3. The 4+4 coefficients are 
then obtained by inverting a simple matrix equation, relating the 
rank-2 tensor components to rank-3 tensor components. 

Owing to the use of periodic boundary conditions in the hor- 
izontal directions, the z-component of the mean magnetic field 
is conserved and equal to the initial value, i.e. B3 — 0. Therefore 
the value of 0133 is here of no interest. 

It is convenient to discuss the results in terms of the quanti- 
ties 



7 = 2"( a 21 - "12) 

Vt = |(?7ii +7722), 

5 = - ?7i2)- 



e 7 = \(a 2 \ + 0-12) 

£r, = |(77ll - 7722), 



(21) 
(22) 
(23) 



Furthermore, the remaining or otherwise important coefficients 
are analyzed individually. The most important of these are the 
diagonal components of and 7721. The former are responsible 
for the generation of magnetic fields in helical turbulence and 
the latter can drive the mean-field shear-current dynamo in non- 
helical turbulence with shear (Rogachevskii & Kleeorin 2003 
120041 

To normalize our results, we use isotropic expressions of a 
and rjt as obtained from first order smoothing, i.e. 



a 



In k' 1 

3 u, rms"'f 1 



(24) 



where the root mean square velocity is a volume average and the 
Strouhal number, 



St 



T c U r 



s kf, 



(25) 



has been assumed to be of the order of unity. In order to actually 
compare our results with those of FOSA, anisotropic expressions 
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Table 1. Summary of the runs. The numbers are given for the statistically saturated state. Here, k = fe/fci, Ma = u Tms /(gd) 1 ^ 2 , 
and L H = L x = L y , 



Run 


grid 


Pr 


Ra 


Rm 


Pm 


Sh 


Co 


e 


Ma 


k 


L H 


A 


128 3 


1.37 


3.1 


10 h 


37 


5 








- 


0.046 


1 


4 


B 


128 3 


1.37 


3.1 


10 5 


35 


5 





0.36 


0° 


0.043 


1 


4 


C 


128 3 


1.37 


3.1 


10 5 


46 


5 


-0.14 





— 


0.058 


1 


4 


D 


128 3 


1.37 


3.1 


10 5 


35 


5 


-0.18 


0.36 


0° 


0.044 


1 


4 


Al 


128 3 


1.37 


3.1 


10 b 


37 


5 








- 


0.046 





4 


A2 


128 3 


1.37 


3.1 


10 5 


37 


5 








- 


0.046 


1 


4 


A3 


128 3 


1.37 


3.1 


10 5 


38 


5 








- 


0.048 


2 


4 


A4 


128 3 


1.37 


3.1 


10 5 


38 


5 








- 


0.048 


3 


4 


Bl 


64 x 128^ 


1.37 


3.1 


10 b 


33 


5 





0.38 


0° 


0.042 


1 


2 


B2 


128 3 


1.37 


3.1 


10 5 


35 


5 





0.36 


0° 


0.043 


1 


4 


B3 


256 x 128 2 


1.37 


3.1 


10 5 


32 


5 





0.40 


0° 


0.040 


1 


8 


B4 


128 3 


0.69 


6.1 


10 b 


1.6 


0.1 





0.32 


0° 


0.049 


1 


4 


B5 


128 3 


0.69 


6.1 


10 5 


3.2 


0.2 





0.32 


0° 


0.050 


1 


4 


B6 


128 3 


0.69 


6.1 


10 5 


7.8 


0.5 





0.33 


0° 


0.049 


1 


4 


B7 


128 3 


0.69 


6.1 


10 5 


16 


1 





0.32 


0° 


0.050 


1 


4 


B8 


128 3 


0.69 


6.1 


10 5 


32 


2 





0.32 


0° 


0.050 


1 


4 


B9 


128 3 


0.69 


6.1 


10 5 


75 


5 





0.34 


0° 


0.047 


1 


4 


BIO 


256 3 


0.69 


6.1 


10 5 


155 


10 





0.33 


0° 


0.049 


1 


4 


Bll 


128 3 


1.37 


3.1 


10 b 


35 


5 





0.36 


0° 


0.044 





4 


B12 


128 3 


1.37 


3.1 


10 5 


35 


5 





0.36 


0° 


0.043 


1 


4 


B13 


128 3 


1.37 


3.1 


10 5 


35 


5 





0.36 


0° 


0.045 


2 


4 


B14 


128 3 


1.37 


3.1 


10 5 


35 


5 





0.36 


0° 


0.045 


3 


4 


B15 


128 3 


1.37 


3.1 


10 b 


35 


5 





0.07 


0° 


0.044 


1 


4 


B16 


128 3 


1.37 


3.1 


10 5 


33 


5 





0.15 


0° 


0.042 


1 


4 


B17 


128 3 


1.37 


3.1 


10 5 


35 


5 





0.36 


0° 


0.044 


1 


4 


B18 


128 3 


1.37 


3.1 


10 5 


33 


5 





0.78 


0° 


0.041 


1 


4 


B19 


128 3 


1.37 


3.1 


10 5 


29 


5 





1.74 


0° 


0.037 


1 


4 


B20 


128 3 


1.37 


3.1 


10 5 


20 


5 





6.43 


0° 


0.025 


1 


4 


B21 


128 3 


1.37 


3.1 


10 b 


35 


5 





0.36 


0° 


0.044 


1 


4 


B22 


128 3 


1.37 


3.1 


10 5 


36 


5 





0.35 


15° 


0.045 


1 


4 


B23 


128 3 


1.37 


3.1 


10 5 


37 


5 





0.34 


30° 


0.046 


1 


4 


B24 


128 3 


1.37 


3.1 


10 5 


38 


5 





0.33 


45° 


0.048 


1 


4 


B25 


128 3 


1.37 


3.1 


10 5 


37 


5 





0.34 


60° 


0.047 


1 


4 


B26 


128 3 


1.37 


3.1 


10 5 


40 


5 





0.32 


75° 


0.050 


1 


4 


B27 


128 3 


1.37 


3.1 


10 5 


43 


5 





0.29 


90° 


0.054 


1 


4 


CI 


128 3 


1.37 


3.1 


10 b 


43 


5 


-0.03 







0.054 


1 


4 


C2 


128 3 


1.37 


3.1 


10 5 


42 


5 


-0.06 







0.052 


1 


4 


C3 


128 3 


1.37 


3.1 


10 5 


46 


5 


-0.14 







0.058 


1 


4 


C4 


128 3 


1.37 


3.1 


10 5 


66 


5 


-0.19 







0.083 


1 


4 


Dl 


128 3 


1.37 


3.1 


10 b 


37 


5 


-0.03 


0.06 


0° 


0.046 


1 


4 


D2 


128 3 


1.37 


3.1 


10 5 


37 


5 


-0.07 


0.15 


0° 


0.043 


1 


4 


D3 


128 3 


1.37 


3.1 


10 5 


37 


5 


-0.18 


0.36 


0° 


0.044 


1 


4 


D4 


128 3 


1.37 


3.1 


10 5 


37 


5 


-0.36 


0.73 


0° 


0.044 


1 


4 


D5 


128 3 


1.37 


3.1 


10 5 


37 


5 


-0.83 


1.66 


0° 


0.038 


1 


4 



need to be used. Such expressions have been computed in the 
past (e.g. Radler 1980; see also Kapyla et al. 2006a) and are 
given for the a-effect, 7, and rj t by 



(0) 

XX 

(0) 

vv 

(0) 



a 
a 

7 

(0) 



-2r c u x d y u 



-r c d z u 2 z , 

T C U„ , 



(26) 
(27) 
(28) 
(29) 



where we have used integration by parts and assumed that r c 
does not depend on spatial coordinates. The correlation time can 
be presented in terms u rms and fcf by assuming a value for St. 



2.3. Averaging and error estimates 

In the present study a mean quantity is considered to be a hori- 
zontal average, defined via 



— 1 



1 T 



F(x + x',y + y',z,t)dx'dy'. (30) 



Except for special terms such as the shear terms in Eqs. and 
(fT3l l. this formulation corresponds to simple horizontal averag- 
ing (for details see Brandenburg et al. 2008a). An additional time 
average over the statistically steady part of each simulations is 
also applied. The fluctuating magnetic fields b v,q are reset to 
zero after periodic time intervals in order to avoid the complica- 
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tions arising from the growth of these fields; see the more thor- 
ough discussions in Sur et al. ( 120081 ) and Mitra et al. ( 120091 ). 

We estimate errors by computing the standard deviation a 
for each depth and dividing this by the square root of the num- 
ber of independent realizations N of the dynamo coefficients. 
We consider the time series between two resets of the field b p ' q 
to represent an independent realization. For a typical run, N is 
between five and ten. 

2.4. Corresponding mean-field models 

In order to determine how well the derived dynamo coefficients 
describe the dynamos seen in direct simulations of Paper I, we 
construct a one-dimensional mean-field model where the test 
field results can be used directly as inputs. We start from the 
mean-field induction equation, Eq. (Q~|), which can be written us- 
ing the vector potential 



Ai = -Uj ti Aj + (XijBj - (rjij + r]6ij)(j,oJj 



(31) 



where the dot on Ai denotes a time derivative and Uj^Aj = 
(SA y , 0, 0) is the shear term. The mean magnetic field is given 
by B = (— A y , A x , 0), the mean current density is given by 

fioJi = —Ai , and primes denote z-derivatives. The coefficients 
a.ij and rjij are taken directly from the test field results leaving 
little freedom in the model. We can, however, turn on and off 
any component of o;^ and rjij when needed in order to study the 
effects of the different coefficients individually. 

3. Results 

In a similar fashion as in Paper I we perform four types of sim- 
ulations which we label as follows: in set A neither rotation nor 
shear is present whereas in set B rotation is added. In set C only 
shear is present, and finally in set D both rotation and shear are 
used. Parameters such as the strengths of rotation and shear, as 
measured by Co and Sh, respectively, are varied within each set 
to probe the parameter space. Summary of the runs is presented 
inTableQ] 

The fluid Reynolds numbers in our simulations are quite 
modest so we cannot consider our flows to be highly turbulent. 
However, the flows are irregular enough to remain time depen- 
dent in all cases, as can also be seen from various animation^ 

3. 1. Set A: no rotation nor shear (Co = Sh = 0) 

The simplest case we can consider with the present setup is one 
with no rotation and no shear. In that case no net helicity gen- 
eration or a-effect are expected. However, due to the density 
stratification, the turbulence is inhomogeneous. This can lead to 
a non-zero pumping, or 7-effect (e.g. Krause & Radler [T980l ). 

The horizontally and temporally averaged transport coeffi- 
cients from Run A with Co = Sh = and Rm « 37 are 
presented in Fig. Q] The results show that the kinetic helicity 
is small and that the mean values of the diagonal elements of 
ctij are of the order of 0.1 ao with errors clearly larger than 
the mean. Vanishing diagonal elements of is in accordance 
with expectations from symmetry arguments. There is however 
a non-zero pumping effect directed upward (downward) in the 
lower (upper) part of the convectively unstable layer. The sign 
of the pumping is inconsistent with the diamagnetic effect, i.e. 
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Fig. 1. The three topmost panels show the time-averaged vertical 
profiles of kinetic helicity, a xx , and a yy , respectively. The fourth 
and fifth panels show 7 with e 7 , and r/t with e n , respectively. The 
lowermost panel shows S (solid line) and r/ yx (dashed). From 
Run A with Co = Sh = 0, and Rm 37. The shaded areas 
between the thinner lines indicate error estimates as described in 
Sect. 12. 31 The vertical lines at z = (0, d) denote the base and top 
of the convectively unstable layer. 
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7 oc —d z u 2 z (e.g. Radler IT9681 I and differs from earlier results 
from convection simulations using the imposed field method 
(Ossendrijver et al. 120021 Kapyla et al. |2006a) and other di- 
agnostics (e.g. Nor dlund et al. [19921 Tobias et al. [19981 1200T1 
Ziegler & Riidiger 2003). However, this result is obtained for 
test fields for which k/ki = 1, whereas the imposed field results 
use a uniform field with kjk\ = 0. For a uniform test field the 
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Fig. 2. Coefficients 7 (top panel) and % (bottom panel) as func- 
tions of k from Runs A1-A4 with Rm « 37, and Co = Sh = 0. 
The dotted lines show the FOSA results for 7 and 77, according 
to Eqs. and (O, respectively, with St w 1.59. 



pumping effect indeed changes sign and is thus consistent with 
the earlier numerical studies and the diamagnetic effect (see the 
upper panel of Fig. The FOSA-prediction, Eq. (|28T > for the 
turbulent pumping is in qualitative agreement with the simula- 
tion result for k/k\ = but opposite to the results for k/k\ 
greater than that. 

At first glance the magnitude of the turbulent diffusivity 
seems quite high: the maximum value is more than six times the 
isotropic reference value rjto suggesting that St rs 6. However, 
the high value of rj t turns out to be related to the normalization: 
if an anisotropic expression, i.e. Eq. ( |29l , is plotted alongside 
rjt the Strouhal number is roughly 1 .6, not six, for our standard 
case k/ki — 1, see the lower panel of Fig. (0. The profile of the 
turbulent diffusivity coincides with that of the vertical velocity 
squared as predicted by Eq. d29i >. When k is increased, the pro- 
file of rjt stays roughly the same and the magnitude diminishes 
roughly in proportion to fc _1 . The quantities e 7 , e n , 5, rj xy , and 
rj yx are compatible with zero in all runs in set A. 

3.2. Set B: only rotation (Co ^ 0, Sh = 0) 

When rotation (corresponding to the north pole, 9 = 0) is added 
to the system, non-zero negative kinetic helicity is produced due 
to the fact that g ■ f2 < 0. Although the a-effect is not directly 
proportional to the helicity in the anisotropic case of stratified 
convection, it can still be a useful proxy. Figure [3] shows the 
results for Run B with Co »s 0.36 and Rm 35. We find that the 
diagonal components of ay are positive in the upper part of the 
convectively unstable region where the helicity is most negative. 
The negative maxima at the base of the convection zone are, 
however, not reflected by the kinetic helicity. The results in Fig. [3] 
were obtained for k/k% = 1. The profiles of a xx and a yy are 
more in line with the profile of the helicity for k/k\ = (see 
Sect. l3.2.3l for more details on the fc-dependence). 

In comparison to Run A, the pumping coefficient 7 shows a 
deeper maximum in the upper half of the convection zone and 
somewhat decreased value in the lower half. The profile and 
magnitude of the turbulent diffusivity are similar to those in the 
nonrotating case. The coefficients rj yx and r\ xy are equal in mag- 
nitude and of opposite sign. This leads to a positive (negative) S 
in the convection zone (overshoot layer) with magnitude peak- 
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Fig. 3. Same as Fig.Q] but for Run B; Co 

and Rm ss 35. 



0.36, 9 = Sh = 0, 



ing close to twice ry t o- The quantities e 7 and e v are small, as 
expected from symmetry arguments 

3.2.1 . Dependence on horizontal system size 

The profiles and magnitudes of the two diagonal components 
of a are very close to each other in our standard case (Run B) 
shown in Fig. [3] Although the number of convection cells in the 
domain is quite small (of the order of ten or less), the isotropy 
of a is a preliminary indication that the horizontal system is still 
large enough to give representative results relevant for a larger 
ensemble. If the system size is too small, the derived turbulent 
transport coefficients may no longer be meaningful (Hughes & 
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Fig. 4. From top to bottom: kinetic helicity, a xx , a yy , 7, %, and Fig. 5. From top to bottom: kinetic helicity, a xx , a vy , 7, %, and 

5 as functions of horizontal system size from Runs B1-B3. The S as functions of Rm from Runs B4-B10. The linestyles are as 

linestyles are as indicated in the lowermost panel. Co « 0.36, indicated in the lowermost panel. Co « 0.33, 9 = Sh = and 

8 = 0, Sh = and Rm w 35 in all runs. Re w 15 in all runs. 



Cattaneo 2008 ). In order to study the convergence of our results, 
we have performed simulations with three box sizes where the 
horizontal extent Ln = L x — L y is either 2c?, 4d, or 8d, re- 
spectively, and the vertical extent of the box is kept unchanged. 
These runs are labeled (from the smallest to the largest) as Bl, 
B2, and B3, where B2 is the same as Run B (cf. Fig. [3]). All three 
runs are relatively slowly rotating with Co « 0.36, Sh = 0, and 
Rm ss 35. The results are shown in Fig. [4] It is obvious that the 
differences between the runs are very small and the two larger 
systems are virtually identical. The only statistically significant 
difference is the anisotropy of a in the convectively unstable 



region for Run Bl with the smallest system size. We can thus 
be fairly confident that the standard box size with — id is 
sufficiently large. This is consistent with the results of Hughes 
& Cattaneo (2008 ) who found that the ratio Ln/d needs to be 
larger than two for the value of a to be reasonably representative 
of the system. 

3.2.2. Dependence on Rm 

One of the basic expected properties of turbulent dynamos is 
that they should be 'fast', i.e. the growth rate of the dynamo, 
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Fig. 6. From top to bottom: a, 7, ry t , and <5 as functions of k for 
runs B11-B14 with Rm m 35 and Co « 0.36. Linestyles as 
indicated in the lowermost panel. 



and thus the transport coefficients, should not depend on the 
molecular magnetic diffusion provided that Rm ^ 1. Figure 
[3] shows the transport coefficients as functions of Rm for fixed 
Re w 15, Co w 0.33, and 9 = Sh = from Runs B4-B10. 
The magnetic Reynolds and Prandtl numbers vary in the ranges 
1.6 . . . 155 and 0.1 ... 10, respectively. We find that the trans- 
port coefficients show no statistically significant dependence on 
Rm for Rm J> 8. The only appreciable departures occur for the 
two lowest Reynolds numbers. However, our definition of the 
Reynolds number depends on k{ (= 2n/d) which is not as well 
defined as in, e.g., forced turbulence simulations, so the coeffi- 
cients may depend on 77 for values somewhat larger than Rm = 1 
in the present case. These results agree with those obtained for 
isotropic turbulence (Sur et al. 2008; Brandenburg et al. 2008a 1 
Mitra etal. 120091 




Fig. 7. From top to bottom: kinetic helicity, a, 7, r/ t , and 5 as 
functions of rotation from Runs B15-B20. The linestyles are as 
indicated in the second panel from the bottom. 9 — 0, Sh = 0, 
and Rm w 20 . . . 35. 



merically (see the two previous sections) and combine the data 
of the a-effect into a single coefficient a. The same applies to 
the FOSA expressions, Eqs. d26b-(l27]i. which we combine into 



y(0) 



(0) 

XX 
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). For k/ki — 0, the a-effect shows a 



3.2.3. Dependence on wavenumber k 

The results for nonrotating convection (see Fig. [U indicate that 
at least the pumping effect can experience not only a change 
in magnitude but also a qualitative change when the wavenum- 
ber of the test field is varied (for corresponding details see 
Brandenburg et al. 2008b ). It is of great interest to study whether 
similar effects can occur for the a-effect. Our results for the stan- 
dard case of Rm w 35 and Co « 0.36 from runs Bl 1-B14 are 
shown in Fig. [2] According to symmetry arguments, the diagonal 
components of ay are the same for 9 = 0. We confirm this nu- nonrotating case. 



more uniform positive value in the convectively unstable region 
than for kjk\ = 1. For larger k the negative region in the deeper 
layers shifts towards the top and the maxima of the profile di- 
minish. The pumping coefficient and turbulent diffusion behave 
very much the same as in the nonrotating case, cf. Fig. [2] The 
coefficient S diminishes rapidly as k increases. The FOSA ex- 
pression for the a-effect is in qualitative accordance with the 
k/ki = result, but fails to capture the details of the profile 
(see also Kapyla et al. 2006a). The Strouhal number required to 
match 77^°' with ij t is 1.37 which is somewhat smaller than in the 
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Fig. 8. From top to bottom: kinetic helicity, a xx , a yy , 7, r/ t , and 5 
as functions of colatitude 9 from Runs B21-B27. The linestyles 
are as indicated in the lowermost panel. Co « 0.29 — 0.36 and 
Rm « 35 — 43 in all runs. 



3.2.4. Dependence on Co 

The a-effect from rotating convection simulations in setups sim- 
ilar to ours has been studied in numerous papers in the past using 
the imposed field method (e.g. Ossendrijver et al. 120011 [2002; 
Kapyla et al. 2006a). Our results for the Coriolis number depen- 
dence of the kinetic helicity, a, 7, r/ t and S from Runs B15-B20 
are shown in Fig. [7] We find that the kinetic helicity and a in- 
crease monotonically as functions of rotation in accordance with 
the results of Ossendrijver et al. (2001). The vertical pumping 



effect shows little dependence on rotation although the Coriolis 
number changes by two orders of magnitude. 

Interestingly, the turbulent diffusivity shows a marked de- 
crease for rapid rotation. The coefficient 5 is positive in the con- 
vection zone and negative in the overshoot layer for slow ro- 
tation. The magnitude increases rapidly until Co w 0.2, after 
which (5 changes sign near the top. This negative region increases 
with rotation. Similar results, i.e. monotonically decreasing rj t 
and a first increasing and then decreasing <5 were obtained from 
forced turbulence simulations by Brandenburg et al. (2008a). 

The combined effect of increasing a and decreasing 7y t sug- 
gests that the large-scale dynamo was possibly subcritical in 
the runs with only rotation in Paper I and other earlier studies 
(e.g. No rdlund et al. 1 19921 Brande nburg et al l 19961 Cattaneo & 
Hughes 120061 Tobias et al. |2008l l, but that it could be excited 
for more rapid rotation. The validity of this conjecture is given 
some credibility by Kapyla et al. (2009a) who find clear large- 
scale dynamo action for Co > 4 for a similar setup as used here 
and in Paper I. More detailed discussion of these results can be 
found in the aforementioned reference. 



3.2.5. Dependence on 

The latitude dependence of the coefficients for Co w 0.29 — 0.36 
and Rm w 35 — 43 from Runs B21-B27 is shown in Fig. [8] The 
colatitude is varied from 0° (north pole) to 90° (equator) in in- 
crements of 15 degrees. The kinetic helicity and the diagonal 
components of <Xy decrease monotonically towards the equator. 
The latter are approximately equal and show a similar latitude 
dependence. This is consistent with the results of Ossendrijver et 
al. ([2002) and Kapyla et al. (2006a) with a comparable Coriolis 
numbefl Vertical pumping also decreases monotonically from 
the pole to the equator. The equatorial profile of 7 is quite simi- 
lar to the nonrotating run, cf. Fig.Q] The trend is similar to that 
seen in earlier studies, e.g. Kapyla et al. (2006a). The variation 
of ?7t is smaller than for the other components, but a weak in- 
creasing trend from the pole towards the equator is seen. The 
coefficient 5 shows a clear decreasing trend as a function of co- 
latitude and is consistent with zero at the equator - in accordance 
with symmetry considerations. 



3.3. Set C: only shear (Co = 0, Sh ^ 0) 
3.3.1. Simulation results 

The next case to consider is that of shear only. We use uniform 
shear of the form U = (0, Sx,0), where S < 0, resulting in 
g ■ (V x U) > 0, and positive kinetic helicity, as expected. See 
Fig. [9] for representative results from Run C with Sh « —0.14 
and Rm « 46. We find that introducing shear into the system 
produces an anisotropic a-effect. The profile of the a xx compo- 
nent is similar, but of opposite sign, to that in the case with only 
rotation; see, e.g., Fig. [3] The magnitude of this component is 
also quite large, i.e. up to twice the estimate o.q already for rather 
weak shear of Sh —0.14. The a yy coefficient, relevant for dy- 
namo excitation, is positive, but the magnitude is only about one 
fifth of a xx . Moreover, the error bars are so large that the value 
is hardly statistically significant. These results demonstrate that 
linking the a-effect to the negative of the kinetic helicity can 
be misleading. Figure [TUl shows the coefficients as functions of 
shear for Runs C1-C4. We have to restrict the study to rather 



" The definition of Coriolis number in the present study is smaller by 
a factor of 2n in comparison to previous studies. 
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Fig. 9. Same as Fig.Q] but for Run C with no rotation and just 
shear; Co = 0, Sh « -0.14, and Rm « 46. 



modest values of Sh because shear, in the absence of rotation, 
promotes generation of large-scale vorticity (e.g. Elperin et al. 
2003; Kapyla' et al. 2009b). Although our largest value of Sh is 
still rather modest, the a xx component is quite large, up to three 
times oto. The a yy component, however, remains small and pos- 
itive for all values of Sh without a consistent trend as a function 
of shear. 

The turbulent pumping in Run C has a similar profile as in 
the cases with Co = Sh = (Run A) and Co ^ (Run B) with 
k/ki = 1 with downward pumping near the surface and upward 
pumping in the lower part of the convectively unstable region. 
The profile and magnitude of the turbulent diffusivity is also very 
similar to previous cases. The pumping effect and turbulent dif- 



3 

> 

3 



0.25 

0.20 

0.15 

0.10 

0.05 

0.00 
-0.05 




-0.5 



0.0 



0.5 



1.0 




o 



0.5 
0.0 
-0.5 
-1.0 

6 

O 

f? 4 

\ 

e- 2 
o 



0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 
-0.6 































V s y / '■ - 









-0.5 



0.0 



0.5 



1.0 



Sh = -0.19 (C4 
Sh = -0.14 (C3 
Sh = -0.06 (C2 
. Sh = -0.03 (CI 










..V 










• N - - - ' 










\ \ 










" - -v.- / 





-0.5 0.0 0.5 

z/d 



1.0 



Fig. 10. From top to bottom: kinetic helicity, a xx , a yy , 7, r/t, and 
r\ yx as functions of shear from Runs C1-C4. The linestyles are 
as indicated in the second panel from the below. Co = and 
Rm w 42 — 46 in all runs. 



fusivity are decreased when the magnitude of Sh is greater than 
0.06. The results for increasing a and decreasing r/ t as functions 
of shear are opposite to those obtained from helically forced tur- 
bulence with shear (Mitra et al. 2009 ). However, the comparison 
for the a-effect should be done with caution because in Mitra et 
al. ( 120091 ) a arises essentially due to the external forcing and is 
only modified by the action of shear whereas in the present case 
a is due to the interaction of shear, stratification, and turbulence 
themselves. 

The rj yx component, which can drive a mean-field shear- 
current dynamo for rj yx S > 0, is of interest because it can pro- 
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Fig. 11. Growth rates A from a one-dimensional mean-field 
model with transport coefficients from Runs C1-C4 presented 
in Fig. [10] Linestyles as indicated in the lower panel. 



vide an explanation for the dynamos seen in recent dynamo sim- 
ulations (Paper I; Hughes & Proctor 2009). In the present case 
where S < 0, r\ yx should be negative to excite the shear-current 
dynamo. There appear to be consistently negative regions of r\ yx 
at the interface of the convectively unstable region and the over- 
shoot layer, and in the upper layers of the convection zone. The 
upper negative region is more pronounced for Sh = —0.06 and 
Sh = —0.14. However, the errors of these quantities are of the 
same order of magnitude as the mean value, cf. the bottom panel 
of Fig. [9] These results tend to agree with earlier findings from 
forced turbulence (Brandenburg et al. 2008aj Mitra et al. [2009 ) 
where r\ yx for the most part was positive or compatible with zero. 

3.3.2. Mean-field dynamo models 

In Paper I clear large-scale dynamo action was found from a sim- 
ulation with Sh w —0.08 whereas for Sh w —0.03 the solution 
was marginal. In the range —0.08 > Sh > —0.22 the growth 
rate of the large-scale field was proportional to the modulus of 
the shear parameter S. Bearing these results in mind and having 
obtained the turbulent transport coefficients for the correspond- 
ing parameter regime we are in a position to apply the test field 
results in a mean-field dynamo model. 

Using the full test field results for and rjij in the corre- 
sponding mean-field model indicates that a dynamo is excited 
for — Sh > 0.14, see the growth rates of the large-scale field 
presented Fig. [TT] It is interesting to study what are the rel- 
ative importances of the different effects: first we turn off the 
off-diagonal components of rjij in which case the magnetic field 
is generated by the a-effect and the shear-current effect is ab- 
sent. We find that the growth rate decreases but is still positive 
for the same cases as before. On the other hand, a "pure" shear- 
current dynamo, i.e. where ctfjj = r\ xy = 0, is also excited for the 
same runs with a very similar A as in the a-shear case. In com- 
parison to the simulations of Paper I, we find that the growth 
rates from the mean-field model are consistently significantly 
smaller. These results and the fact that no dynamo was found 
for — Sh = 0.06 would seem to indicate that an incoherent a- 
shear dynamo (e.g. Vishniac & Brandenburg |199Tb is also op- 
erating in the full simulations. However, we should remain cau- 
tious when comparing the direct simulations and the mean-field 
model because the transport coefficients were determined for a 
single value of k whereas many other wavenumbers are available 
in the simulations. 
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Fig. 12. Same as Fig. [JJ but for Run D with both rotation and 
shear; Co w 0.36, Sh rj -0.18, and Rm Rj 37. 



3.4. Set D: rotation and shear (Co ^ 0, Sh ^ 0) 
3.4.1. Simulation results 

When rotation is added to the system where a large-scale shear 
is already imposed, the vorticity generation is suppressed (e.g. 
Yousef et al. 2008b ; Paper I) and it is possible to study higher 
values of Sh. Representative results for Run D with Co w 0.36, 
Sh w —0.18, and Rm Rj 35 are shown in Fig.[T_2] The results 
for the kinetic helicity and a yy seem to behave additively when 
comparing with the runs with only rotation (Run B; Fig.O and 
only shear (Run C; Fig. [9). The a xx component is somewhat 
smaller than in Run C with shear only which is consistent with 
oppositely signed contributions due to shear and rotation. 



12 



Kapyla et al.: Alpha effect and turbulent diffusion from convection 



3 

> 
3 



0.0 
-0.5 
-1.0 



0.5 
0.0 
-0.5 
-1.0 
-1.5 
-2.0 
-2.5 









" \ 










*>7 - - / / 

v-„ y ■ 










\ / 






-0.5 


0.0 


0.5 1.0 










N.J 








<•<■■ ••■ 1 










V-' 






-0.5 


0.0 


0.5 1.0 



Q 




6 

O 

2 




2 

O 
X 

P 



-1 





-0.5 


0.0 


0.5 


1.0 






7 


\ 








-0.5 


0.0 


0.5 


1.0 




... Co=1.66 
-. Co = 0.73 
- - Co = 0.36 
Co = 0.15 
Co = 0.07 


D5) : 

D3 : A?- 
Dl) :/£— 


"•*•>> 














j 




-0.5 


0.0 

z/d 


0.5 


1.0 



Fig. 13. From top to bottom: kinetic helicity, a xx , a yy , 7, i] t , 
and j] yx as functions of Coriolis number for Runs D1-D5. The 
linestyles are as indicated in the lowermost panel. Sh = — i Co 
and Rm « 37 in all runs. 



We note that in a recent paper, Hughes & Proctor (2009) 
found that the a-effect is virtually unchanged when shear is 
added to a rotating system. In their case the shear profile is 
proportional to cosy. The resulting large-scale vorticity is then 

W z oc siny which leads to a\Y^ °c (G • W)Sij oc siny (e.g. 
Radler & Stepanov 120061 ) where G symbolically denotes the 
inhomogeneity of the turbulence. However, Hughes & Proctor 
(2009 ) show a volume average of a over the full upper half of 
the domain in which case the contribution of a^ w ^ cancels out. 
This explains the absence of any modifications of a due to shear 
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Fig. 14. Same as Fig.QTjbut f° r Runs D1-D5 shown in Fig. Qj] 



in their case, but for us this is not the case because for our shear 
profile W z = S — const. 

The profile of 7 is quite similar to the rotating case, i.e. 
Run B; see Fig. [3] with the exception that the off-diagonal com- 
ponents of ctij show considerable anisotropy as manifested by 
the parameter e 7 . The turbulent diffusion shows a profile com- 
mon to all the other simulations, but here the diagonal compo- 
nents of r)ij show evidence of mild anisotropy with e v peaking 
near the middle of the convectively unstable region with a max- 
imum value of e v ks r)to- The quantity S is compatible with zero 
whereas r\ yx exhibits a similar profile and magnitude as 5 does 
in the rotating case, cf. Fig. [3] which indicates that r\ xy w r\ yx , 
as opposed to r\ xy w — r\ yx in Run B. 

The kinetic helicity and the turbulent transport coefficients 
as functions of Co, keeping the ratio — S/Sl = 1 constant, are 
shown in Fig. [13] The helicity is increasing in a similar fash- 
ion as, albeit slower than, in the absence of shear, compare with 
Fig. [7] The components of the a-effect are highly anisotropic 
with the main contribution of a xx being due to shear and that 
of a yy due to rotation (compare with Figs. 171 and [Toll. The value 
of a xx is somewhat decreased in comparison to the cases with 
shear only whereas a yy is almost unaffected. This is in qualita- 
tive agreement with adding the contributions of runs from Sets B 
and C with corresponding Co and Sh, respectively. 

The profiles of 7 and r]t are similar to those with only rota- 
tion. The differences of 7 and rj t as a function of Co are small 
and for the most part fall within the error bars. The profile of 
r) yx is very similar to 5 in the case of only rotation with negative 
values near the base and top of the convectively unstable region 
with positive values in between. The profile and magnitude of 
r) yx remains essentially fixed for Co il 0.15. We note that in the 
simulations with shear and rotation the J~i x J-effect may also 
contribute to the generation of large-scale magnetic fields (e.g. 
Radler [19691 Radler et al. l2003l Pipin et al. l2008l It is, how- 
ever, not altogether clear how to disentangle the transport coef- 
ficients responsible for the shear-current and f2 x J-dynamos in 
the present case. Thus the off-diagonal components of Tfo con- 
tain contributions from both effects in Set D. 



3.4.2. Mean-field dynamo models 



We follow here the same procedure as in Sect. 13.3.21 to study 
dynamo excitation for the Runs D1-D5. Using the full and 
rjij tensors in the one-dimensional mean-field model indicates 
that all of the runs in Fig. [13] are capable of driving a dynamo. 
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Neglecting the off-diagonal components of r)ij decreases A by 
approximately a third whereas for a pure shear-current dynamo 
the growth rate is roughly half of the model where the full Oy 
and r]ij tensors were used. For an a-shear dynamo A is positive 
for all runs except the slowest rotation case with Co = 0.07 
whereas for the combined shear-current and f2 x J-dynamo also 
Co = 0.15 is mildly subcritical. It is interesting to note that the 
simulation results of Paper I fall roughly on top of the uppermost 
line, i.e. where the full a^ and r/ij are used, in Fig. [14] This is in 
contrast to the case of only shear where the growth rates from the 
mean-field model were clearly smaller than in the corresponding 
simulations of Paper I. However, the mean-field model does not 
reproduce the declining growth rate for Co J> 1.1 that is ob- 
served in the direct simulations. 

4. Conclusions 

We obtain turbulent transport coefficients governing the evo- 
lution of large-scale magnetic fields from turbulent convection 
simulations with the test field method. We study the system size 
and magnetic Reynolds number dependences of the coefficients. 
This is important because spurious results can be expected for 
small Reynolds numbers or when the aspect ratio of the do- 
main is too small (Hughes & Cattaneo 2008). We find that for 
our standard system size, L^/d = 4, the coefficients are essen- 
tially identical to those obtained with a horizontal extent that is 
twice as large. As a function of Rm, all the coefficients are es- 
sentially constant for Rm 8. This is in accordance with the 
theory but at odds with results from certain imposed field cal- 
culations (e.g. Cattaneo & Hughes 2006). In these calculations 
the magnetic field is allowed to evolve until saturation which can 
cause strong quenching even if the imposed field itself is weak. 
This is particularly important if closed boundary conditions for 
the magnetic field are imposed, in which case magnetic helicity 
conservation can lead to catastrophic quenching (Vainshtein & 
Cattaneo 1992). More reliable results for the kinematic a-effect 
with the imposed field method can be obtained by resetting the 
magnetic field before it grows too large or develops substantial 
gradients (e.g. Ossendrijver at al. 120021 Kapyla et al. [2006a). 
More detailed comparison of the imposed field and test field 
methods is important, but beyond the scope of the present study. 

The earlier determinations of transport coefficients from 
convection simulations have used the imposed field method (e.g. 
Ossendrijver et al. 1200 II [20021 Kapyla et al. I2006at which yields 
the components of ay but does not deliver r/ij because the im- 
posed field is uniform. The test field method does not suffer from 
this restriction and r/ij and the fc-dependence of the coefficients 
can be extracted. We find that for fc/fci = 0, i.e. for a uniform 
field, the results for a and 7 are consistent with those obtained 
from imposed field calculations, provided the magnetic field is 
reset before it grows too large and substantial gradients develop. 
As k is increased, however, the qualitative behaviour of the co- 
efficients changes. This is indicated by a partial sign change of 
a and a complete sign change of 7; see Figs.[2]and[6] 

The turbulent diffusivity shows a robust behaviour regardless 
of the parameters of the simulations: the profile is proportional 
to the vertical velocity squared, u\, as predicted by FOSA (e.g. 
Radler [T980l l. The value of rj t decreases almost proportional to 
fc _1 , and shows a declining trend as a function of rotation and 
shear. 

For the present parameters, the a-effect increases monotoni- 
cally as rotation is increased. As a function of latitude, the diag- 
onal components of ay have a similar magnitude and peak near 
the pole with declining values towards the equator. The a-effect 



induced by shear is highly anisotropic: the a xx component has 
a similar profile and magnitude, but opposite sign, as a xx and 
a yv in the case of only rotation. This component also increases 
monotonically as a function of shear, whereas the shear-induced 
a yy remains small regardless of the strength of the shear. In the 
runs where rotation and shear are present, the diagonal compo- 
nents of ay are roughly the sums of the corresponding coeffi- 
cients in the cases with rotation and shear alone. 

In addition to the a-effect, the rj yx component can contribute 
to a shear-current dynamo when rj yx S > 0. In our case, where 
S < 0, such dynamo action is possible if r\ yx < 0. We find that 
this coefficient shows negative regions near the base and near the 
top of the convectively unstable region, but the errors are of the 
same order of magnitude as the negative mean values in most 
cases. 

In order to connect to earlier work, we use the test field re- 
sults in a one-dimensional mean-field model in order to under- 
stand the excitation of dynamos using identical setups as in di- 
rect simulations (Paper I). We study here only the cases with 
shear and consider large-scale dynamos in the rigidly rotating 
case elsewhere (Kapyla et al. [2009a). The presently used dy- 
namo model ignores fc-dependence and is therefore likely to 
be too simple to fully describe the large-scale fields in the di- 
rect simulations. Nevertheless, the present results, taken at face 
value, seem to indicate that in the case with shear alone the de- 
rived dynamo coefficients are not sufficient to explain the dy- 
namo but that an additional incoherent a-shear dynamo might 
be needed. This conjecture is based on the fact that mean-field 
a-shear and shear-current dynamos are both excited with simi- 
lar growth rates which, however, are significantly smaller than 
those obtained from direct simulations in Paper I. Furthermore, 
a large-scale dynamo was marginal for Sh = —0.03 in Paper I, 
whereas for Sh = —0.06 it was found to be slightly subcritical in 
the present study. On the other hand, for the case with both shear 
and rotation, no additional incoherent effects seem to be needed. 
We find that in this case the regular a-shear dynamo produces 
larger growth rates than the combined shear-current and f2 x J 
dynamo but neither effect alone seems to be strong enough to 
explain the dynamos in Paper I. 

On a more general level, mean-field dynamo models of the 
Sun and other stars rely on parameterisations of turbulent trans- 
port coefficients. Even today, the majority of solar dynamo mod- 
els bypass this problem and ignore most of the turbulent effects 
and rely on phenomenological descriptions of the a-effect and 
turbulent diffusion that are not without problems theoretically. 
On the other hand, some attempts have been made to incorpo- 
rate the results for the transport coefficients from imposed field 
studies in mean-field models of the solar magnetism (e.g. Kapyla 
et al. |2006bt Guerrero & de Gouveia Dal Pino [20081 ) and mod- 
els employing more general turbulence models have recently ap- 
peared (e.g. Pipin & Seehafer 2009 ). We feel that this is a worthy 
cause to follow further with the present results. 
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